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Abstract 

We construct explicit multivortex solutions for the complex sine-Gordon 
equation (the Lund-Regge model) in two Euclidean dimensions. Unlike the 
previously found (coaxial) multivortices, the new solutions comprise n single 
vortices placed at arbitrary positions (but confined within a finite part of 
the plane.) All multivortices, including the single vortex, have an infinite 
number of parameters. We also show that, in contrast to the coaxial complex 
sine-Gordon multivortices, the axially-symmetric solutions of the Ginzburg- 
Landau model (the stationary Gross-Pitaevskii equation) do not belong to a 
broader family of noncoaxial multivortex configurations. 
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I. INTRODUCTION 



Topological solutions of nonlinear partial differential equations on the plane have been a 
subject of intensive investigations in recent years, with applications ranging from nonlinear 
optics to cosmic strings. The simplest type of topological solitons arise in systems involving 
one complex field. In this case the solitons realise a S 1 —>■ S l map of the circle of a sufficiently 
large radius on the (x, y)-plane into a unit circle in the internal space (Reip, Imip) and can 
be classified according to the Brouwer degree of the map, or the winding number, 

Q = lim -!- i d (axgip) = — / e nk d n d k (aigip) d 2 x. (1) 
R-^oo Zn J Cr Zn J R 2 

One celebrated model of this kind is the Ginzburg-Landau model (also known as the station- 
ary Gross-Pitaevskii equation) arising in the description of boson condensates (in particular, 
superconductors and superfluids) [0,0: 

VV + ^(1 - l^l 2 ) = 0. (2) 

(Here V = id x +jd y .) Another system of a similar type which received considerable attention 
in literature, is the Heisenberg ferromagnet with easy-plane anisotropy ||. In terms of the 
stereographically projected field this can be written as 

In physics literature the winding number (|l|) is usually referred to as vorticity and the 
planar topological solitons as vortices. Both the Ginzburg-Landau and the easy-plane ferro- 
magnet equation are well known to possess axially-symmetric solutions describing n vortices 
sitting on top of each other, ip( n > = e m6l $ n (r), where the function $ n (r) — > 1 as r — ► oo. 
(Here r and 6 are the polar coordinates on the (x, y) -plane.) Of fundamental importance, 
both physically and mathematically, is the question of whether noncoaxial multivortices 
exist. However, despite some encouraging recent insights H], very little is known about this. 
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The present paper is devoted to another physically meaningful systemP] which shares a 
lot of similarities with (fj) and (|3]), the so-called complex sine-Gordon equation: 

The equation (|j) is similar to @ and (§) in that it is also an equation for one complex 
field on the plane, and that it also possesses coaxial multivortex solutions. An important 
difference, however, is that this model is integrable and hence the complex sine-Gordon 
multivortices are given by explicit analytic formulas whereas their Ginzburg-Landau and 
magnetic counterparts are available only numerically, even in the single-vortex case. The 
single- vortex solution of eq.(^) has the form 

^ {1 \r,6) = (5) 

where Iq and ii are the modified Bessel functions. Using a purely algebraic recursive pro- 
cedure || one can construct axially-symmetric solutions of arbitrarily high vorticity. For 
example, the Q = 2 and Q = 3 solutions are given by 

wr - w) 

and 

, (3) (r B) = [h(r) ~ h(r)][I 2 (r) - lf(r)] + ^r)^) - h{r)f 3W 
V {:) [I (r)-I 2 (r)}[Io(r)I 2 (r)-2P 1 (r)+I$(r)} ' [) 

respectively. 

The primary objective of this paper is to show that the complex sine-Gordon equation 
has an infinite-parameter family of exact, explicit, noncoaxial multivortices. This family 
includes the axially-symmetric n-vortex solution which, therefore, admits a continuous into 
n spatially separated single vortices. We calculate the energy (i.e. the Euclidean action) 



1 Some references to the applications of the complex sine-Gordon theories can be found in the 
concluding section below. 
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of the noncoaxial n-vortex configuration and demonstrate that it does not depend on the 
positions of the individual vortices. This implies that the individual vortices making up the 
n-vortex configuration are non-interacting. 

Having established the existence of a general n-vortex solution for equation (^), a natural 
question to ask is: Do the Ginzburg-Landau model (|2]) and the easy-plane ferromagnet ([3]) 
have multivortex solutions other than the axially-symmetric ones? Confining ourselves to the 
case of the Ginzburg-Landau equation, we attempt to answer this question via the analysis 
of the spectra of linearised excitations of its symmetric solutions. If the axially-symmetric 
solution admitted a p-parameter nonsymmetric continuation, the corresponding linearised 
operator would have p zero eigenvalues. We study the linearised spectra of the Ginzburg- 
Landau axisymmetric solutions numerically; the upshot of this study is that they have only 
three zero modes related to obvious symmetries of the equation and therefore do not belong 
to a broader family of noncoaxial multivortices. 

On the contrary, each of the coaxial multivortices of the complex sine-Gordon theory 
admits an infinite number of zero- frequency excitations. This is an expected fact, of course, 
given the existence of the infinite-parameter nonsymmetric configurations. However, as we 
show below, the presence of the infinite number of zero modes in this case can be established 
even without knowledge of the nonsymmetric generalisations. In other words, the possibility 
of the nonsymmetric continuation of the coaxial multivortex could have been predicted 
simply on the basis of the analysis of its linearisation. 

The paper is organised as follows. In the next section we derive the noncoaxial n-vortex 
solution and in section 3 discuss some of its general properties. Section 4 deals with the 
simplest special case of the new solution when it depends only on one free parameter. We 
consider, in detail, the n — 1, n — 2 and n = 3-solutions, and then extrapolate our conclu- 
sions to the situation of the general n. The energies of the multivortices are evaluated in 
section 5. In section 6 we study, numerically, the linearised excitations of axially-symmetric 
multivortices of the Ginzburg-Landau equation and compare them to their complex sine- 
Gordon counterparts. Finally, some concluding remarks are made in section 7. 
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II. THE GENERAL MULTIVORTEX SOLUTION 

As in ||, we start with rewriting the second-order equation (^) as a system of two first 
order equations: 

^ (n - 1} + (1 - W n ~ 1] | 2 ) = 0, (8) 

Q^(n) _ _ |^(")| 2 ) = 0. (9) 

Here d = d/dz, d = d/dz, and z = (x + iy)/2, z = (x — iy)/2. This first-order system 
has a field-theoretic interpretation of its own; it is nothing but the Euclidean version of the 
massive Thirring model ||[7|. Both ip^ 1 *) and ip^ satisfy equation (§]), hence eqs.(||)-(H) can 
be seen as the Backhand transformations relating two solutions of the complex sine-Gordon 
equation. 

Let n — 1 in eqs.([8])-(|9]). For any ip^ eq.([8D is solved by ip^ = 1. Letting = 1 in 
eq.(g), we get 

= 1- |^ (1) | 2 . (10) 
Decomposing ip^ = f + ig, the imaginary part of fllPP yields 

d x g = d y f, 

whence we can define the potential T{x, y) such that / = d x T and g = d y J- ', or, equivalently, 

v> (1) = aF. (ii) 

The real part of fliCf ) is then 

V 2 ^+ (V^) 2 = 1. 
Letting T = lnZ , this reduces to the (modified) Helmholtz equation 

V 2 Z -Z = 0. (12) 
In polar coordinates, the general solution of (|i~2"D, regular at the origin, is given by 



Z (r,9) = J2 I m(r)a m (9), (13) 

m=0 

where I m (r) is the modified Bessel's function of order m (of the first kind); a m = (3 m cos(m#+ 
5 m ), and /3 m ,5 m are arbitrary real constants. (In particular, all (3 m with m greater than a 
certain M can be set equal to zero in which case the series ( |i~3"|) becomes a finite sum.) 
Returning to the variable ^ > we obtain 

7 /.(i) _ r ie f E I'mi^mjO) , j E I m (r)d g a m (9) 1 , v 

where I' m = dl m /dr and all sums run over m = 0,1,. ..oo. We will assume that (3q ^ 
0; otherwise the above solution is singular. Without loss of generality we can let <To = 
flo cos 8q — 1 in ( |lj|) and this convention will be implied throughout this paper. At infinity, 
the solution flTJ]) tends to e ld ; more precisely 

^) (riJ ) = e«{l_i^5 + ei+«!l + ^} as^oc, (15) 

where «, /^i and i/i are functions of 6 1 : 

2<9<?E^n ^(4m 2 -l)(x m 
« = , tn = ^ , vi = -deni{9). (16) 

Therefore, equation (|j"4|) gives a solution with infinite number of parameters and vorticity 



Q = 1. For purposes of this paper we will be referring to (14) as the "general one- vortex 
solution". Now the general solutions with Q = 2 and all higher vorticities can be obtained 
in a purely algorithmic way. We simply use eq.(§) to express ip^ n > via ijj^ 1 ^: 

^ (n) -- 1 _|^(n-i)i2 ^ (n " 1) ' - = 2,3,-.. (17) 
This recursive procedure can be made rather efficient by introducing auxiliary variables 

oo 

Z k (r,9) = J2^mtk+m + l m Zk-m), k = 0, ±1, ±2, .., (18) 



m=0 



where 



Ur,9)=I a (r)e M , lm = %e^. (19) 



Note that Z- k = Z k \ also note that the function Zq has already been defined before (see 
eg . (|l~3"D . ) What makes the variables Z k useful is that the differential operators d and d act 
on them simply as index lowering and raising operators: 

dZ k = Z k -u dZ k = Z k+1 . (20) 

Indeed, writing d as e~ l9 (d r — -dg), d as e t9 (d r + -dg) and using the identities || 

d-^rra Im—1 Im+X Tnlm Im—1 Im+1 



dr 

one can readily verify that 



(21) 



1)g f(Urn m \ i{m-l)6 _ t , 

y dr r / 



fdl m m \ _ i(m+i)e _ t 

—, J-m — J-m+l e — im+1- 

\ dr r / 



From here the relations (|20f) are straightforward. 

Recalling that T = \tlZq and using equation (|TTD , the general 1-vortex solution (|i~4"D can 
be written simply as 

^\r,9) = ^. (22) 



Now applying ( |TTD and using ( |20"D gives 



^ 2 Hr,e)= *}-? 22 2 , (23) 



^,0)^ _ ^3(^0 - Z-^Zt) + £-i£ 2 2 + g| - 2Z a ZiZ 2 ^ 



and so on. Setting to zero all /3 m with m > 1, eqs.(|22]), ( |23| ) and (|24"D reduce to the axially- 
symmetric vortex solutions (|), (§) and (0). 

The explicit expressions of the multivortices with n > 4 become cumbersome and we 
restrict ourselves to producing only their asymptotic behaviours as r — ► 00: 
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,<».(, e)e - = 1 ^± + eM + , (JL K(e) + ^) + o (-L ) , (25) 



where /x n and z/ n are defined by recurrence relations 

4/!„_i 

A*n — fJ-n-2 H T H — ) 

n — 1 2 

n+1 ds/^-i m 

= T^n-l : 71 - 2 ' ( 26 ) 

n — 1 n — 1 4 

with /i = and /ii, z/i and k as in (|i~6l) . Equations (|25|)-(^6|) can be easily proved by 
induction with the help of the Backlund transformation The recurrence relation for 

/i n can be easily resolved yielding an explicit expression 

/j n = n 2 /xi + — L doK, n>2. (27) 
(Unfortunately, there are no similar closed formulas for v n .) 



The relations ( 18[) can be seen as expansions over the eigenfunctions of the angular 



momentum, with f3 m being the coefficient of the eigenfunction associated with the orbital 



quantum number m. Accordingly, solutions (|22|) -(|25l) can be interpreted as orbital defor- 



mations of the axially-symmetric multivortices. Below, in section IV, we will discuss several 
particular orbital deformations in more detail. 

III. SOME GENERAL PROPERTIES 
A. Regularity and convergence considerations 

It is not difficult to realise that inequality 

oo 

^ 1 ( 28 ) 

m=l 

is sufficient to ensure the regularity of the general 1- vortex solution ( p2|) . (It is not necessary 
though; see the next subsection.) Indeed, due to the positivity of I m (r) for r > and 
the second identity in (PT|), we have Io(r) > hmij) and I\{r) > hm+iij) for all m > 1. 
Combining these inequalities with Io(r) > I\(r) gives /o( r ) > Im{ r ) for all m > 1 and 
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r > 0. This latter inequality, taken together with (|28|), guarantees that Zq > and hence 
the solution ( p2|) is bounded on the entire (x, y)-plane. 

In case of infinitely many nonzero coefficients j3 m we need to make sure that the series 



in (28) converges. This can be accomplished by imposing, for example, that 



lAnl < q r 



(29) 



with some < q < 1. The inequalities fl28[) and (|2~9"D are also sufficient for the convergence 
of the series in the asymptotic formula (|i~5|). 

Next, the 2- vortex solution resulting from the Backlund transformation (|17|) will only be 
regular if the seed 1-vortex solution is bounded by 1 in absolute value, that is, if \Zi\ < Z . 
That the latter inequality holds true can be easily verified using the representation 

i r as 

' -G(£)e z£+ ^\ U{Z) = 

m=0 



Zt 



2m J £ k+l 
\e\=i 



(30) 



which arises by replacing £ m in eq.([L^) by 

U{r,e)=I m {r)e vmd 



2- f j^e Xp [zl + 1 ). 



(31) 



Eq. fl3"T|) , in turn, follows from the integral formula for the modified Bessel function ||: 

= C^ eXP {l( C + f)}' 

ICI=i 

where we only need to set ( = £e 10 . (We remind the reader that z = \re l9 .) To show 



that \Zi\ < Z , assume that the series ^2 \Pm\ converges, with eq.(|28|) being in place. The 
function G(e lLp ) = 1 + J2m=i @m cos(rmp + 5 m ) is obviously positive for all ip and hence using 
equation (|30|) with £ = e itp we get, finally, 



\Zi\ = — 
1 11 2vr 



2tt 



d(pe^ v G(e %v )e 



1 



2- 



< — / dcpG(e l *)e rcos ^ +G 
2ti 



Zn. 



Q.E.D. 
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B. Translations of the vortices 



It is interesting to note that for certain choices of parameters /?i , /?2 5 - - - solutions ip^\ 
ip( 2 \ etc., describe pure translations of the corresponding coaxial multivortices. The 
translations are associated with infinite numbers of nonzero /3's. Consider, for example, the 
translation along the x-axis: 

x = x — R, y — y. (32) 



The transformation (|32|) can be written as one formula, 

f cos(# + <f) = r cos(6 l + if) — R cos <p, 

which holds for an arbitrary fixed angle (p. Hence, picking up G(e llfi ) = e~ RcoS(fi in the 
representation (p0|), we obtain 

Z k = — [ d(pe- ikv e- RcoBV e ram ^ + ^ = &(f, 6), (33) 
2tt J 
o 

where ^ were defined in eq. (|19]) : £k(r,9) = Ik{r)e lke . According to (|33|), in the reference 
frame (x, y) our functions Zk of which the solution is to be built, have the form flT8| ) with 
all f3 m = except (3$ = 1. In other words, in the translated reference frame the solutions 
ip( n ) constructed using G(e ltp ) = e~ Rcosip , have the form of coaxial n- vortices. 

The orbital coefficients (3 m associated with the translation can be found by expanding 
the function G = e~ RcoSip in the Fourier series: 



oo 

- R cos ip 



J2 hme^ + 1^) , (34) 

m=0 

where 

7o = ^o(i?); 7m = (-l) m / m (i?), m>l. (35) 

It is not difficult to check that the series ([33]) converges. Indeed, the large-m asymptotic 
behaviour of the modified Bessel's function is given by Horn's formula [§]: 
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I m (R) = -^=expjm (l + In - (m + In mj (l + O Qj) > m ^ oo. (36) 

Using (|36|), one can easily verify that the coefficients (|35D pass the ratio test: 

lim = 0. (37) 

We can normalise the coefficients according to our convention that (3q be equal to 1. This 
is done simply by replacing 7 m in eq.(|35|) with 



7o = ^; 7 m = (-l) m ^, m>l. (38) 

(Note that the corresponding /3 m , (3 m = 2(— l) m I m (R)/ I (R), do not satisfy the sufficient 
condition (|28|) . Despite that, the translated vortex is perfectly regular.) 

Thus we conclude that the recursion procedure (|17D-(|18D with an infinite sequence of 
nonzero orbital coefficients 7 m defined by eq.(^), gives rise to the coaxial multivortex centred 
at the point x = R, y = 0. In the next section we will show that choosing a finite number 
of nonzero 7's may also result in a shift of the vortex; however that shift will always be 
accompanied by a deformation. On the contrary, the infinite sequence (^) produces a pure 
translation. 

Our final remark in this section is on the convergence of yet another series: 

00 

m=0 

(This is a useful by-product of eq.(|3T|).) The fact that the above series converges for all 
k allows one to use the asymptotic formula (|15|) in the case of the translated 1-vortex. 
Restricting ourselves to terms of order r _1 , we find that the axially-symmetric (undeformed) 
vortex (|]) moved to the point (x — R, y — 0), has the asymptotic behaviour 

^ M)= e«{l-l + ^ + (I)} „™ (39) 
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IV. ONE-PARAMETER DEFORMATIONS OF THE AXI ALLY- SYMMETRIC 

MULTIVORTICES 



A. The n = 1 vortex 

In this subsection we analyse in detail the simplest situation of a single nonzero orbital 
perturbation: /3& = (3 7^ for k equal some fixed m, and (3% = for all other k. Without 
loss of generality we can set 5 m = and consider (3 to be non-negative: < (3 < 1. Eq.(p2|) 
gives 

(1) = /i(r)e i0 + (g/2) [J m+ i(r)e^ m + 1 ) e + / m _ 1 (r)e^( m ~ 1 ) e ] 
^ ' ~ I (r)+pi m (r) C os(m9) ' [ ] 

The first three orbital perturbations (m = 1,2,3) of the one- vortex solution are shown in 
figure 0. The unperturbed vortex (i.e. eq.(|40D with f3 = 0) is also reproduced for comparison. 

The one-parameter solution ( flUD is symmetric with respect to the rotation 9 — > 9 + 2-K/m 
in the (x, ?/)-plane. This accounts for the number of symmetric folds in the modulus of t/^ 1 -* 
seen in figures |l](b-d). The solution (fHJ) with m > 3 is different from figures |l](b-d) only in 
that it will have m symmetric folds. 

It also follows from the cyclic symmetry that out of all one-parameter perturbations, 
only the m = 1 perturbations give rise to the shift of the vortex from the origin. This 
shift breaks the rotational Z m -symmetry completely and therefore, is compatible only with 
m — 1. This can be easily verified by the Taylor's expansion at the origin, 

V> (1) (r,0) = ye"' 51 +0(r). (41) 

According to (0), the value of ip^\ r _ is not equal to zero — unless (3i = 0. 
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FIG. 1. (a): the modulus squared j?/^ 1 )! 2 of the unperturbed n = 1 vortex, (b-d): its m = 1, m = 2 
and m = 3 orbital deformations. Note that in panel (b) we used eq.(^C|) with 8\ = tt/2 for the better 
visualisation. This gives rise to the y-shift of the vortex from the origin (and not the x-shift as in the 
case of <5i = 0.) Similarly, in panel (d) we set 63 = n. 



B. The n = 2 multivortex 

The one-parameter deformations of the coaxial two- vortex with m = 1,2,3 are shown 
in figure 0. Each of these three perturbations splits the repeated zero of tp^(x,y) into 
several single zeros, or, equivalently, splits the coaxial two- vortex into several monovortices. 
The perturbations with m = 1 and m = 2 give rise to two single vortices. For m = 2, 
the indices of the newly born zeros are equal due to the discrete rotation symmetry 9 — > 
9 + 2n/m = 9 + 7r; hence in this case the perturbed solution consists of two individual 
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vortices of vorticity Q\ = Q2 = +1. Here the vorticity of the k-ih vortex is defined as the 
index of the corresponding zero of the field: 



Q k = ^l d(argV0, (42) 




FIG. 2. The modulus squared of the n = 2-solutions. (a) the coaxial 2-vortex (/3 = 0); (b) the 
m = 1 perturbation (here 5i = 7r/2); (c) m = 2; (d) m = 3. In (b,c,d) we set (3 = 1 to ensure the 
maximum possible separation of individual vortices. 
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FIG. 3. The level curves of |^^ 2 ^| 2 for the coaxial two-vortex (a) and its first three orbital perturba- 
tions: m = 1 (here 5% = it/ 2) (b); m = 2 (c); m = 3 (d). In each of these plots = 1. Only level 
curves with sufficiently low values of \i/j\ 2 are shown for visual clarity. 

For m = 3 there are four zeros. They are not clearly visible in figure [|(d) but become 
evident if one plots the level curves of |V^| 2 > figure [| First of all, we have an antivortex 
with Qq = — 1 sitting at the origin. In addition, there are three vortices with vorticities 
Qi,2,3 = +1 placed symmetrically around it. The fact that the origin is a zero with index 
— 1, can be readily concluded from the Taylor expansion for small r: 



V> (2) (r, 



4 



-^e~ i9 + 0{r 2 ). (43) 
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Perturbations with higher orbital numbers do not produce the splitting of the 2-vortex. 
For all m > 4 the solution has a zero only at the origin. For example, for m = 4 and 5 the 
Taylor expansions about r = are: 

~ 2 ' (3 



# 2 >(r,i 



m=4 



- 1 - '-e 



+ 0(r 4 



(44) 



,2i6 



+ 0{r s 



(45) 



m=5 



Since < (5 < 1, the argument of the term in brackets in ( [44]) is a periodic function of 9 
and hence the index at the origin equals +2. The sole effect of perturbations with higher 
orbital numbers on the 2-vortex solution amounts to the symmetric deformations, similarly 
to the effect of m > 2-perturbations on the 1-vortex solution. 



C. The n = 3 multivortex 

We conclude our discussions of one-parameter deformations with the case of the (n = 3)- 
solution. The absolute values of ip^ are shown in figure [| for m = 1, 2, 3 and in figure |6| 
for m = 4 and m = 5. (As we explain below, all higher orbital numbers cannot produce the 
splitting of the 3-vortex and are less interesting, therefore.) The location of the positions 
of individual vortices is facilitated by plotting the level curves of l?/^ 3 -*! 2 , figures [5] and [7]. 
A poor visibility of zeros in figure |6|(b) is due to a significant disproportion in the widths 
of the central and surrounding vortices. In figure 0(b) we have attempted to improve the 
visualisation by contouring the central vortex at lower levels of l^l 2 than the surrounding 
ones. 
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FIG. 4. The modulus squared of the (n = 3)-solutions. (a) the coaxial 3-vortex {j3 = 0); (b) its 
m = 1 deformation (here S% = — vr/2); (c) m = 2; (d) m = 3. In (b-d), (3 = 1. 
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FIG. 5. Level curves of \tp\ 2 for the coaxial three-vortex (a) and its lowest orbital deformations, (b) 
m = 1 (here ft = -vr/2); (c) m = 2; (d) m = 3. In (b-d) (3 = 1. 



For m > 2, the indices of the surrounding zeros can be found from the symmetry consid- 
erations provided the index of the zero at the origin is known. The Taylor expansion about 
r = gives: 



^ {3) (r, 



m=2 



^ {3) (r, 



= ^e-« + 0(r 3 ) 

771=4 



^ (3) (r,#) 



m=5 



/3r 2 
~16 



e - 2 ^ + 0(r 3 ). 



(46) 
(47) 
(48) 



For higher orbital quantum numbers, m > 6, the (n = 3)-solution consists of just one 
vortex (of vorticity +3) sitting at the origin. 
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FIG. 6. The modulus squared of the m = 4 (a) and m = 5 (b) orbital deformations of the coaxial 
three- vortex. Here j3 = 1. 
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FIG. 7. The level curves of |^| 2 for the two solutions shown in figure ||. (a) m = 4; (b) m = 5. 
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D. The one-parameter deformations of the coaxial multivortices with general n 



Having analysed a number of particular combinations of n and m, we can now formulate 
general conjectures on the indices of one-parameter multivortex solutions. We have identified 
three distinct cases depending on the relation between n and m. 

(i.) The simplest situation occurs if m > 2n. In this case, it is sufficient to determine 
the index of the vortex at the origin. We conjecture that the Taylor's expansion, as r — >• 0, 
is 



4 



(«) 



2 n n\ 



1 - (-l) n S m on^e- 2ine 



And i /r\( n-\-\ 



e 



0{r n+1 ), m > 2n. (49) 



(Here 5 mi 2n is Kronecker's delta.) According to (f49|) , in this case we have a single multivortex 
sitting at the origin, with vorticity Qq = +n. 

(ii.) Next, let n < m < 2n. In this case we suggest the following expansion, as r — > 0: 

R(—~\ \ n+l r .m-n 

^ (n) = ^T^T7 M e " (m ~ n)0 + O(r n+1 ), n<m<2n. (50) 

2 m ~ n+l {m — njl 

Eq.(pOl) shows that there is a vortex with vorticity Qq = — (m — n) at the origin and m 
vortices with indices Qi,2,...m = +1 placed symmetrically about the origin in accordance 
with the rotation symmetry 9 — » + 2-n jm. 

(iii.) Finally, it remains to consider the case m < n. This is the most nontrivial situation 
as in this case there may be no vortices at the origin at all. [See e.g. figures H(b),(c) and 
U(b),(d)]. The analysis of several combinations of n and m suggests that the index Qq at 
the origin is the smallest-modulus remainder (positive or negative) from the division of n 
by m. That is, Qq is equal to the smallest of all positive and negative integers q satisfying 

n = km + q, \q\ < m, (51) 

where A; is a positive integer. In case there is a positive and negative remainder of equal 
modulus (i.e. if m = 21 and n = km±l with / > 0), the index equals the positive remainder: 
Qq = I. Each of the km surrounding vortices has index +1. These are placed symmetrically 
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about the origin (except for the case of m = 1, of course) in accordance with the rotation 
symmetry 9 — > 9 + 2n/m. For instance, the n = 4, m = 2 solution consists of four vortices 
centred on the x-axis and grouped in two pairs symmetrically with respect to the origin. 

All particular solutions considered so far verify equations (f4l^-(|51~|). However we do not 
yet have proofs of these formulae in the case of general n and m. 



V. THE ENERGY OF THE MULTIVORTICES 

In this section we calculate the action integral for the vortices. In the phenomenological 
theory of phase transitions this integral characterises the free energy of the system. (Ap- 
parently for this reason this quantity is also referred to as energy in mathematics literature 
|i~0Hi~2| — although it gives the stationary Hamiltonian only for particular (2+l)-dimensional 
extensions of our (2+0)-dimensional model (p.) Besides its fundamental role in physical 
applications, the action can be used as a powerful variational tool in a purely mathematical 
analysis of stationary and moving topological solitons (see e.g. ||10|-[T3|]). 

The action integral for equation (f|) has the form 

E W = / I W + 1 - W) = / m±™* + 1 - A, (52) 



Dj.yl-W 2 Jd r \ 2(1- 



where the integration is over a disc of a large radius R centered at the origin. The massive 
Thirring model (B)-(H), in its turn, is derivable from the action 



^[^ (n - i) ,^ (n) ] = / iv^- 1 )^^ - v (n) <M n_i) + (i^ (n_i) r - i)(^ (n) r - 1) + c.c.jd 2 x. 

J D R 

(53) 

The two systems are equivalent; expressing t^™ -1 ) through from eq.(^) and substituting 
into (^) produces eq.(^) for tp( n \ Accordingly, the integrands in ( |52"D and (^) differ only by 
a divergence: 

5[^ (n - 1} ,^ (n) ] = 2E[iP {n) ] -I V ■ A {n) d 2 x; n > 1, (54) 

JDt, 



where 
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A (n) = -Vln(l - |^ n) | 2 ) + 2W(|^ (n) | 2 ) V x arg^ 



(n) 



(55) 



and 



1-p 



+ ln(l -p). 



Alternatively, we can express through x ); substituting this into (§) produces eq.( 
for f/^™ -1 ). This equivalence is, again, reflected by the corresponding actions. We have 



where 



S'[V' (n-1) ,V (n) ] = 2£[V> (n ~ 1) ] - / V • B^- 1 )^; n>2, 

Jd r 



B (n-i) = _vin(l - |V> (n_1) | 2 ) - 2W(|^ (n ~ 1) | 2 ) V x arg^^ 1 ). 



(56) 



(57) 



Both A^ n ) and B^™ -1 - 1 are regular in the finite part of the (x, y)-plane and therefore the 
double integrals in fl54|) and (|56|) can be transformed to contour integrals over the boundary 
of the disc Dr. 



^[^ {n_1) ,^ (n) ] = 2E[^ {n) ] - & A (n) -rdl, n > 1, 

JdD R 



(58) 



S'[V' (n_1) ) V' (n) ] = 2E[^ n ' 1] ] 



Here f = r/r. For large r eq.([25|) gives 



B^-rdZ, n>2. 



|^(n)|2 = l_5+ 1 



.-, ; 1 + k 2 \ n(n — 1 

n 2 2/ii + — - — ] + — -Kg 



O 



A J 2 

As one can easily check using flTB]) , the 0(r~ 2 ) term in (^) is a total derivative: 

l + K 2 



+ 



-^ 2 ln£ 



(J,, 



and so (^) can be written as 



|^)| 2 = l-- + i,%AK 



(59) 



(60) 



(61) 



Therefore 
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W(|^ (n) | 2 ) = --ln- + — 1 + 1-) asr^oo, (62) 

n n n \r J 

and using this in (|58| ) - ([59|) we obtain, finally, 

S[^ n - l \^ n) ] = 2E[^ n) ] -4tt ^R-nln^ -n + ij , n > 1, (63) 

^[^ (n_1 \^ (n) ] = 2^[^ {n ~ 1) ] +4vr ^i?- (n - l)ln— ^ -n + , n > 2. (64) 

Subtracting ( |64]) from (^) we arrive at the formula relating the sine-Gordon actions of 
the solutions with vorticities n and (n — 1): 

£(«)-£(«- 1 ) = 27r(2#-nln--(n-l)ln — 2n+l); n > 2. (65) 

\ n n — 1 / 

This equation can be used to calculate, recursively, the energies of all vortices; the only 
outstanding ingredient is the action of the vortex with n = 1. The E^ 1 ' can be found from 
(|63|), provided we know the Thirring action S(ip^ \-ip^). This action is obtained directly 
from eq.(^) where we only need to set ip^ = 1: 

5(^ (0) ,^ (1) )= / + di; w )d 2 x = I + ^e~ ie ) dl. (66) 

Jd r JdD R ^ ' 

Making use of the asymptotic expansion fl25|), eq.(^) yields 



S(^°\ ^ (1) ) = 2tx(2R - 1) + 0{R~ 1 ). (67) 

Finally, the sine-Gordon action of the single-vortex solution is 

= 27r(2 J R-lni2- 1) + 0(R- 1 ). (68) 

The formulas (^) + (^) provide the actions for vortices with any n. As one could have 
expected, the actions do not depend on any of the parameters /3 m ,5 m of the solution (|T4"D. 
Note also that we can use equations (^) and (^) to define the energy of the vorticity- 
free state ip^ = 1. Consistently with one's physical intuition, eq.(|B"5]) with n — > 1 yields 
= 0. 
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Solving the recursion relation ( |65|) with initial condition fl68|) we can obtain the action 
_£:("■) j n closed form: 

n-l 



E {n) = 2tt 



2nR-n 2 lnR + n{lnn - 1) + 2 fc(lnfc - 1) 



k=l 



n>2. 



(69) 



Note that eq.(|69"D remains valid for n = 1 (in which case one should simply disregard the 
sum Y^kZi) an d for n = (in which case one should also set nlnn = 0). 
For completeness, we also evaluate the Thirring action. Using ( |69"D in 

n-l 

2(2n - - 2n(n - 1) In R + 4 fc(ln fc — 1) — 1 



(n-l) ^,(n) 



2tt 



fc=i 



we have 
n > 2. (70) 



Like the previous formula, this expression remains valid for n = 1. 

The main conclusion of this section is that the action (or "energy") of an n- vortex so- 
lution does not depend on parameters /3 m , 5 m and therefore, on the relative positions of 
the individual vortices. This implies that in any (2 + l)-dimensional (i.e. time- dependent) 
extension of the planar complex sine-Gordon theory, vortices may form non-interacting con- 
figurations. 

VI. COMPLEX SINE-GORDON VS GINZBURG-LANDAU: 
ANALYSIS OF ZERO MODES 



The aim of this section is to find out whether the Ginzburg-Landau axially-symmetric 
vortices admit non-symmetric deformations similar to those arising in the complex sine- 
Gordon equation. Let ip^ be an axially-symmetric solution of equation (0): ?/A n )(r, #) = 
$ n (r)e . If this i/j( n > is a member of a broader family of solutions parameterised by p 



continuous parameters ax,a 2 , 



§ n (r)e ine = ip(r,9; a u . . . , a p ) 



(71) 



tti=oc2=...=a p =0 

the equation (^) linearised about the symmetric multivortex will have p solutions of the form 



<%M) 



da. 



j = l,...,p. 



(72) 



ai=a2=—=a p =0 
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Therefore our strategy will be to examine the spectrum of linearised perturbations of the 
Ginzburg-Landau vortices. We will also be considering the linearisation of the complex sine- 
Gordon equation (^); the comparison of the corresponding sets of zero modes for the two 
systems will give rise to some interesting observations. 



A. Linearised eigenvalue problem 

Technically, it is convenient to treat the linearised boundary-value problem as an eigen- 
value problem. With this purpose in mind, we consider the (2 + 0)-dimensional Ginzburg- 
Landau equation (||) as a time-independent reduction of a (2 + l)-dimensional Higgs-field 
equation 

^tt - v 2 v - (1 - IVf)V> = 0. (73) 

(Alternatively, we could have considered it as a reduction of the Gross-Pitaevskii equation 
iipt + V 2 ^ + (1 — IV'I 2 )^ — 0, but the relativistic generalisation ([T3|) is computationally 
advantageous as it gives rise to a symmetric eigenvalue problem.) In a similar way we can 
define a (2 + l)-dimensional generalisation of the planar complex sine-Gordon (|j): 

^ " V ^ - 7^777^ - ^(1 - iVf) = 0. (74) 



(Note that (|74|) is not relativistically invariant and can hardly claim any physical relevance; 
we are introducing this equation just for auxiliary purposes here.) Assuming a solution of 
the form 

i/>(r,6,t) = $„(r)e infl + 6if>{r,6,t) = [$ n (r) + e0(r, 9) cos ut] e ind (75) 

and linearising (|73| ) in small e, we obtain 

1 



V 2 - ^(d 6 + m) 2 - + 2$^0 + ^ = u 2 <P, (76) 



where V 2 = 4> rr + r 1 r . In a similar way, the linearisation of equation ( [74]) gives 
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1 2 2$' n $ 2m 1 



+ 



v 2 A 



1 _$ 



2^2 



+ 2$; 



+ 



r> 1 



1_$ 



2^2 



<j) = id (j), 



(77) 



where the prime over <£> n denotes the derivative with respect to r. Equations ( f7E|) and 
(f?7|) can be regarded as eigenvalue problems, with u; 2 being an eigenvalue and (0, 0) the 
associated eigenvector. Expanding in the Fourier series in 9: 



<j>{r,9) 
and transforming to 



J2 <Pm(r)e im6 = W m (r)+*b m (r)}e ime , 



(78) 



m=— oo 



m=— oo 



«mW = a m + a_ 



(79) 



we obtain a sequence of one-dimensional eigenvalue problems, one for each value of the 
azimuthal number m: 



n,m 



2 / 



(80) 



The operator £„ )m is defined by 
r = 



V2 I n 2 +m 2 



3K(r) - 1 



2mn 



Iran 



-V 2 + ^ + $2(r) - 1 



(81) 



in the Ginzburg-Landau case, and by 
C 



-V 2 + 8„(r)£ + ^ + £„(r) + S) n (r) 



■"n,m 



mMn(r) 



m2l n (r) 

L -L 

' dr r 2 



-V 2 + »„(r) £ + ^ + <t n (r) - S n (r) 



(82) 



in the complex sine-Gordon case. In (|8~2"D we have introduced the notations 

2n 1 



Mr 



r 2 1 — (f)2 

n 

2n $ 2 

7> X — <J) 2 



2$„$ 



n^n— 1 j 



£ n (r) = 2$ 2 - - 1 + 



l + $, 2 , 2n$„_i^ 



+ 



r 2 1 — $ 2 r 1 — $ 2 



r 1 - $ 



2 ' 



(83) 
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The imaginary parts of the Fourier coefficients m satisfy the same eigenvalue problem (|80l 
with £ njm as in (RT1) or (|32"D, where the eigenfunctions (« m ,t> m ) should only be defined by 



u m (r) =b m - 6_ m , v m (r) = b m + b_ m . (84) 

(Note that the eigenfunctions (u m , v m ) do depend on the vorticity n but we are omitting the 
corresponding subscript to keep the notation simpler.) 

Thus the spectrum of linearised excitations of the symmetric multivortex ipW is given by 
eigenvalues of the operator £„ )m with m varying from — oo to oo. In fact in view of equations 
( |79p and flgip we can restrict ourselves to nonnegative m only If (u m ,v m ) is an eigenvector 
of the operator C n>m associated with an eigenvalue u 2 , then (u m , —v m ) is an eigenvector of 
the operator C n _ m associated with the same eigenvalue ui 2 . 

We solved the eigenvalue problem (j80|)-(j8l]) numerically. Before discussing the results of 
the numerical analysis, it is instructive to compare the spectrum structure of the Ginzburg- 
Landau operator fl81]) with that of its complex sine-Gordon counterpart. 



B. The spectrum structure 



For any positive n and m > we introduce two bases of solutions of the linear system 
flSOp-flgip. One basis can be defined by the asymptotic behaviours at the origin. To this 
end, we rewrite the system in terms of u — (u m + v m )/2 and v = (u m — v m )/2 (where the 
subscript m is omitted for simplicity of notation): 
/ 



V 



L m+n + 2$^ 

$2 



$2 

n 



L r 



+ 2$: 



\ 






( ~\ 




u 


= u 2 


u 


/ 


V) 




v) 



(85) 



Here L s = —V 2 + (s 2 /r 2 ) — 1. The advantage of the formulation (|85|) is in that for r — > 0, the 
cross-coupling potentials $ 2 ~ r 2n are small and the equation for u(r) decouples from the 
equation for v(r). Expanding each of the u(r) and v(r) in power series in r (supplemented 
by logarithmic terms where necessary), e.g. 

u{r) = (u + u In r)r p + {ui + Ui In r)r p+2 + ... 
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and substituting into fl85|), one can easily verify that for each pair (n,m) there are four 
solutions with the following asymptotic behaviour: 



^ r m+n [l + o(r)] ^ 



( 



Z-2 



( Q(^ r \m-ri\+2n+2^ \ 



J 



\ 



„|m— n\ 



[1 + o(r)} 



^ r -(m+n)[l + ( r )] ^ 



V 



/ 



o(r~ 



|m— n|+2n- 



-\m— n\ 



[1 +o(r) 



(for m ^ n); Z4 



J 



0(r 2n+2 lnr 
In r ■ [1 + o(r)] 



^-(m+n)+2n+l^ 

(for m = n). 



(Here Z denotes the column (u, v) T , of course.) Transforming back to u m and v m and 
introducing the notation Z = (u m , v m ) T , we have four linearly independent solutions of the 
system ©-(0) (as r -> 0): 



Zj =r ™+«[i + ( r )] 



Z 4 = r 



r l"»-"l[l+o(r)] 



Z 3 = r-( m+ ^[l + o(r)] 



-|m— n| 



[1 + o(r)] 



-1 



(for m ytz n) 



(for m = n) . 



(86) 



Z 4 = In r ■ [1 + o(r)] 

v~v 

The solutions Zi, Z2 are bounded and Z3, Z 4 unbounded for all m and n. 

In a similar way one can show that the linearised complex sine-Gordon fl80|)- (|82D also 
has two bounded and two unbounded solutions near the origin. In this case one should only 
take into account that as r — > 0, 



„, 2n 
21= — + 



2n 



„2n-2 



r 2 (2 n n\y 



+ 0(r zn ) 



2n 



r (2 n n\) 2 



r 2n-l + Q^. 



2n+lN 



n 



(t=—~l + 



2n 2 



„2n-2 



0(r 2n ), 2 = 0(r 



2n\ 



v (2 n n\) 2 

The four basis solutions are given by the same equations fl8"6| ) as in the case of the Higgs 
field. We should emphasise here that eqs.(^) are valid for all u (including u = 0), both in 
the Higgs and sine-Gordon cases. 
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The second basis is denned by the asymptotic behaviours as r — > oo. Consider, first, the 
Higgs system (|S0])-(|8~T|) and let < u < y/2. The four solutions are given by 

/ 



Y, 



3,1 



\ 2iu)r \r z I I 

i ± 2ra2 +L 1/4) - m2 + o ( N 

+ 



2^ 



(87) 



(88) 



The solutions Yi, F 2 , and Y4 are bounded and Y" 3 unbounded as r — > oo. Similarly, for 
all < oo < 2 the linearised complex sine-Gordon fl80D -(j82|) has three bounded and one 
unbounded solution as r — > oo: 

/ 



Yi, 



V 



l-2m 2 
16 



Y 3 4 = e ±v/CTr r p 



/ 

/ 2m 2 -l-p(p+l) | /n / 1\ \ 



V 



2r 



(89) 



(90) 



/ 



2 

r 


1 

4 r 3 


n 

" 2^ 


+ o 


2 

r 


1 


+ o 






2n 


1 


3n 


2 - 


r 


2T 2 " " 


4r 3 




2n 


1 


3n 


2 - 


r 


2T 2 " " 


4r 3 



where p = =p2n(4 — u; 2 ) -1 / 2 — 1. In the derivation of (jS9|)-(|90|) we made use of the asymptotic 
expansions, as r — > oo, of the coefficient functions in (|83|): 

2l„(r) 
®„(r) 

0„ 1 Q„ / 1 \ 

(91) 

Each of the solutions Z^ (i — 1, ...4) can be expanded over the basis Y^: 

i 



For all < w < \/2 in the Higgs case and all < lu < 2 in the complex sine-Gordon case, the 
linear combination T^' m ^Zi(r) — T^' m ' Z 2 (r) represents a solution which is bounded both 
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as r — > and r — > oo. Therefore, in both cases small nonzero uj belong to the continuous 
spectrum. The question that is of concern to us, of course, is whether uj = belongs to the 
continuum; in other words, is there a bounded solution for uj = 0? Surprisingly, the answers 
for the Higgs (alias Ginzburg-Landau) and complex sine-Gordon linearisation, are different. 

Let us start with the Higgs case and let r — > oo. Two asymptotic solutions, Y^, are 
given by eq. (|88f) where we only need to set uj = 0. The other two asymptotic formulas, 
eq. fl87D , cannot be used for uj = 0. Instead, we have two solutions with the asymptotics 

^ m n i f 1 \ 



V 



l + 0(^r) 



Y 1 = hxr-[l + 0(l t )] 



V) 



n = [i + o (M 







V 



(form ^ 0); 



(form = 0). 



(92) 



Therefore, for uj = we only have two, not three, solutions bounded as r — > oo: Y" 2 
and Y 4 . The solutions Z x and Z 2 , bounded at the origin, can still be expanded over the 
basis Yj (j = 1,...4); however, this time in order for the linear combination C\Z\ + c^Zi 
to remain bounded as r — > 00, the constants c\ and c 2 have to satisfy two conditions: 
ciTii' m \o) + C2T2" (0) = and ciT^' m \o) + c^T^ (0) = 0. This imposes a requirement 
on the matrix T^ n ' m \uj) at the point u; = 0: 

T^ m) (0) T 2 ( 3 n ' m) (0) 

which is not a priori satisfied for all values of n and m. Consequently, the zero mode (i.e. 
a solution u m (r), v m {r) pertaining to uj = and bounded for all < r < 00) can only arise 
for some pairs (n,m). 

Consider now the linearised complex sine-Gordon (|80| ) + (|82| ) . As in the Higgs case, the 
asymptotic solutions Y^^r) (eqs.(|9~0"D) are valid for uj = whereas the formulas ( j39|) are not. 
Instead, for uj = the two solutions Yi^if) have the asymptotics 



0. 



(93) 



r 



Y 2 



rn _i_ /rt I In r \ 

2r W Kl^J 

1 + m 2 n^ + C» ) 



/ 



OO. 



(94) 
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Thus in the complex sine-Gordon case we have three bounded solutions as r — > oo: Y±, I2, 
and Y4 (and not two as in the Higgs case.) Accordingly, there always is a solution bounded 
for all r G [0, 00): (u m , v m ) = T^' m) Z 1 (r)-T^ m) Z 2 (r), and so we may conclude that there 
is a zero mode for any pair (n, m) — without having to find this solution analytically or 
numerically. In this sense, in the complex sine-Gordon case the value uo = belongs to the 
continuous spectrum whereas in the Ginzburg-Landau case zero modes can only appear as 
discrete eigenvalues. In the latter case we had to resort to the help of computer. 



C. Ginzburg-Landau zero modes: a numerical search 

To find eigenvalues and eigenfunctions numerically, we replaced derivatives in eq . (jSlD 
(and eqs.(|9^)-((5(J) below) with the second order-accurate finite differences. Typically we took 
a grid with 1000 points on an interval < r < 80 (i.e. with the step size Ar = 0.08). The 
coaxial multivortices $ n (r) of which the potentials in fl5I|) are formed, were pre-computed 
using Newton's method. 

The appropriate boundary conditions at the origin follow from the asymptotic behaviours 

u r (0) = v r (0) = 0, for m ^ n ± 1; 
u(0) =v(0) = 0, for m = nil. (95) 

At infinity, the boundary conditions are chosen to accomodate the bounded solutions in 

u r — > 0, v r — > as r — ► 00. (96) 

Numerically, zero modes appear as small nonzero eigenvalues and one still has to distin- 
guish them from small eigenvalues arising from the continuous spectrum when the infinite 
line is approximated by a finite interval [0,80]. The genuine zero modes can be discerned 
by considering the large-r behaviours of the associated eigenvectors. The "true" zero modes 
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are allowed to decay exponentially (as in fl88|)) or as r~ m (as in (|92|)) whereas the continuous 
spectrum solutions will generically decay as l/\^r (see eq.fl87|)). 

Before calculating the eigenvalues of the Higgs vortices, we tested our numerical scheme 
on the complex sine-Gordon operator (0). As expected, we obtained one zero mode for each 
n and m. (We tested m = 0, 1, 5 for each of n = 1, 2, 3.) The associated eigenfunctions 
were found to coincide with the derivatives of the general n-vortex solution ip^ w.r.t. the 
azimuthal parameters /3 m (see the next subsection.) 

Proceeding to the Higgs system, we examined the one-, two and three- vortex solutions, 
i.e., n = 1,2 and 3. The azimuthal quantum number of the analysed perturbation ranged 
from m = to 5 in each case. Zero eigenvalues were only found for m = and 1; both result 
from obvious symmetry properties of the Ginzburg-Landau/Higgs equation. 

For m = 0, the complex eigenfunction 0(r) of the operator (|76|) associated with the 
numerically- found zero eigenvalue, was found to coincide with the function i$„(r) (i.e., 
Mo = 0, Vq — 2$ n .) This zero mode is related to the U(l) invariance of the Ginzburg- 
Landau equation (g): 



oa 



a=0 



(97) 



For m = l, the numerical eigenvector («i,i>i) associated with the zero eigenvalue was 
found to be equal to (<3>' n , — ^$ n ). (There is, of course, a zero mode for m = — 1 as well: 
(u_ 1 ,v_i) = ($' n , ^$„).) The corresponding complex perturbations (|75|) are the translation 
modes, one <\>e luQ given by d x [$ n (r)e m6 '] and the other one by d y [<3> n (r)e m61 ] . 

Both the translational and the [/(l)-zero modes are well known to workers in this field. 
(See e.g. |Il,|9|jn| ) . There is a simple analytical argument showing that zero modes cannot 



arise for m > 2n pT| ; this fact is also known to specialists. However, the nonexistence of 
zero modes for 1 < m < 2n does not seem to have appeared in literature before. 

Since the only zero-frequency excitations of the axially-symmetric vortices are those as- 
sociated with the phase shifts and translations, we conclude that the coaxial multivortices 
of the stationary Ginzburg-Landau model do not admit continuous nonsymmetric defor- 
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mations. In particular, two vortices sitting, symmetrically, on top of each other cannot be 
continuously separated. This does not mean, of course, that the Ginzburg-Landau model 
does not admit noncoaxial multivortex solutions at all. Multivortex configurations with 
finite separations between individual vortices may exist (and in fact there are indications 
that they do exist [Q). However, the intervortex separations will only admit discrete sets of 
values, or will be allowed to vary continuously but be bounded from below by certain finite 
distances. 

This is in sharp contrast with the coaxial multivortices of the complex sine-Gordon 
theory which can be continuously split and moved apart. As we have shown in previous 
sections, each axially-symmetric solution of this model is a member of an infinite-parameter 
family of solutions corresponding to a specific choice of azimuthal deformation parameters: 
Pi = Pi — ■ ■ ■ — 0. We conclude this section by demonstrating that each of these continuous 
parameters gives rise to a zero mode in the spectrum of the corresponding linearised operator 



D. Zero modes of the complex sine-Gordon multivortices 

The derivative of the general n-vortex solution w.r.t. the azimuthal parameter (3 m (m > 
1), is given by 



dP 



/3i=/3 2 =...=0 1 k=1 _ n ° £ h 



(98) 
ft=/3 2 = ...=0 



To calculate the derivative dip^ /dZk, we notice that each of the multivortices can be written 
as a rational function, 

* w = £55. 0*0 

where the numerator and denominator are homogeneous polynomials in Zx- n , 2 2 - n , ■■■,Z n , 
of degree P (P < 2 n_1 ). These polynomials have the form 
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■A/"^ ^ ^ Ci 1 i 2 ...ip ^i\+i2+...+ip,n Z{^ Z{ 2 . . .Z{ p , 

1— n<ii,i2,—ip<n 

M {n) = 



Di 1 i 2 ...ip 3ii+i2+...+ip,0 Z% 2 ■ --Zip > 



(99b) 
(99c) 



l—n<ii,i2,...ip<n—l 

where C^^...^ and D ili2 ^ ip are real coefficients and is Kronecker's delta. Equations 
( j99|) are straightforward from the recurrence relation fll?]) . Note that in eq. (|99b|) , products 
Z il Z i2 ...Z ip have their indices summing up to n; we will be referring to this property by 
saying that the polynomial J\f^ has level n. In this sense, the polynomial Ai^ (eq.( |99cj )) 
has level 0. When a polynomial of level I is differentiated w.r.t. Zk, its level is lowered down 
to I — k. Therefore, the derivative dip^ / dZk is a rational function whose numerator is a 
polynomial of level n — k and denominator is a polynomial of level 0. Since Z g \a =0 = 
£s = I s { r )e lse , one can easily check that 



dV> (n) 



dZk 



/3 1= /3 2 = ...=0 

with some real function 0^ (r). Substituting into ( |9~8|) we obtain the zero-frequency eigen- 
functions of the operator (|77|) pertaining to the orbital quantum numbers m and — m, re- 
spectively: 



n 



n 



<PU r ) = Yl ^k\ r ) I k+mir)] 0_ m (r) = Y 0fc" ) ('")4-m(r) 



(100) 



fc=l— n fc=l— n 

(m > 1). The zero modes ( |100| ) translate into zero-frequency eigenvectors of the operator 



/ \ 

\v m i 



^ h+ m {r) + h~ m {r) ^ 



k=l—n 



' fc+m 



r) - Ik- m (r) 



(101) 



VII. CONCLUDING REMARKS 

In this paper we have constructed families of explicit solutions of the complex sine- 
Gordon equation. In literature, the complex sine-Gordon model is usually considered in the 
(l+l)-dimensional Minkowski space where it has the form 
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This equation was introduced in the late 1970s in three independent field-theoretic contexts: 
(i) as a reduction of the 0(4) nonlinear cr-model [ Hfl ; (ii) in the description of relativistic 



strings in a uniform antisymmetric tensor field |0J, and (iii) in the theory of massless 
fermions with a scalar contact interaction Later on, eq. (|102|) reappeared in an entirely 



unrelated physical context — it turned out to be equivalent to the Maxwell-Bloch and self- 
induced transparency equations as well as the system governing stimulated Raman scattering 
in nonlinear optics [17|. In mathematics literature it is common to call eq.flTDD the Lund- 



Regge model and write it as 



a xx - ata - Sm o Q (xl ~ Xt) + sin a cos a = 0, 
cos a 



(Xztan a) x = (xttan a) t , (103) 



where sin a and x are the modulus and argument of the complex function ip(x,t) in ( |102| ): 
ip = sina;e* x . Geometrically, the Lund-Regge model gives the Gauss-Codazzi equations for 
the embedding of pseudospherical surfaces into a flat three-dimensional Euclidean space Jl8 . 

More recent studies of the Lorentzian complex sine-Gordon theory ( |102| ) and its Eu- 
clidean counterpart given by eq.(|), were motivated by the fact that eqs. (|102|) and (§) define 
integrable deformations of 2D conformal field theories, more specifically of the SU(2)/U(1) 



coset model and Z n parafermions Other (not unrelated) sources of interest have been 



the search for exactly solvable conformal theories with black-hole background metrics [20 



and exact factorisable ^-matrices on the quantum level ||21|| . In the current mathematics 
literature, the hyperbolic Lund-Regge equation ( |103| ) is being discussed in connection with 
the localised induction hierarchy describing the motion of vortex filaments in an inviscid 
incompressible fluid |2"2|] . Its applications to pseudospherical surfaces continue to attract 
attention (see e.g. pi) while the elliptic equation (H) has been derived in the description of 



the "middle surfaces" of generalised Weingarten surfaces |24 



The (l+l)-dimensional system ( |102| ) was shown to be integrable via the inverse scattering 
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transform [14,25.26] and, by virtue of this approach, broad classes of its exact solutions have 



become available p?| , p8| ,P,[7[] . The inverse scattering framework for the (2+0)-dimensional 
solitons was formulated in ||. Solutions obtained in this way describe nonlinear superposi- 
tions of one-dimensional fronts intersecting at arbitrary angles on the (x, |/)-plane 

In Ref. axially symmetric solutions were constructed describing n vortices (n > 1) 
sitting on top of each other. In the present paper we have shown that each of these coax- 
ial multivortices belongs to an infinite-parameter family of nonsymmetric solutions, which 
includes, in particular, configurations of n single vortices located at separate points of the 
plane. (The constellation of n separated monovortices is not the only splitting possibility 
for the coaxial n-vortex though; there can also be various combinations of vortices and an- 
tivortices with indices summing up to n, the total topological charge characterising that 
particular family.) Our current construction, as well as the one of Ref. 0, is based on the 
Backlund transformation of the complex sine-Gordon equation. (In the axially-symmetric 
case there is also an alternative technique employing the Painleve reduction 0.) For the 
understanding of the structure of the phase space, however, it would be instructive to embed 
the vortex solutions in the inverse scattering formalism; we are planning to return to this 
issue in future publications. The fact that even a one-vortex solution is characterised by an 
infinite number of parameters seems to indicate that vortices should be associated with the 
continuous spectrum rather than discrete eigenvalues (which account for front-like solitons.) 

We have computed the energy (more precisely, the Euclidean action) for each family 
of multivortices. The action is found to be entirely determined by the total topological 
charge and independent of parameters specifying a particular solution within each topological 
class. This is not an unexpected result; the nontrivial dependence of the action on some 
of the continuous parameters entering the solution would contradict the stationarity of this 
configuration. 

That axially-symmetric multivortices of the complex sine-Gordon theory admit noncoax- 
ial continuations, follows already from their spectra of linearised perturbations. As we have 
shown in this work, the existence of a zero mode with any azimuthal number m = 0, 1, 2, 3... 
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is ensured by the availability of the necessary number of bounded asymptotic solutions of the 
linearised equations as r — > and r — > oo. On the contrary, the axially symmetric solutions 
of the Ginzburg-Landau model do not have enough bounded linearised perturbations as 
r — > oo; consequently, the existence of a zero mode for a particular choice of n and m cannot 
be guranteed a priori. Our numerical analysis has demonstrated that the Ginzburg-Landau 
coaxial multivortices have only translational and rotational zero modes and therefore do not 
admit continuous deformations and/or splitting into separate vortices. 
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